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updates 


ALGOL SUBCOMMITTEE FORMED 


On August 22, 1960 the Ancor Working Group held a mecting in Milwaukee. Those attending came from manu- 

fayturers, universities, and computcr-using laboratories. All those attending have written, are writing, or plan to 

ovale Anueot-like compilers or are actively engaged in writing programs in the Ancor 60 language. The attendees, 

“representing 22 organizations, agreed to form a subcommittee of the ACM Computer Languages Committee for the 

_ purpose of maintaining and interpreting the ALGox 60 language. This group is expandable and il is hoped that a Euro- ; 

| pean counterpart may be formed so that actions agreed upon by both groups may be regarded as official interpreta- 
‘tions and changes to Ancor, 60. 

In addition to establishing voting rules (one vote per organization) and electing a chairman (J. Wegstein), the group 

lacted on proposals for changes to ALGou 60 by Peter Naur. There was also a valuable exchange of information on 

- the progress of various compiler implementations, and several papers on implementation techniques were presented. 


bo 


ALGORITHM 16 begin 
CROUT WITH PIVOTING Afi, k] = Ali, k]} - INNERPRODUCT(A[i,p], Alp, kl, 
p, 1, k-1); : 


‘Grorau K. Forsyrus 
ponent Ee Ene err if abs(Afi, kl) > TEMP then 
Atanford University, Stanford, California 3: begin 

TEMP := abs(Ali, k]); imax := i 
real procedure INNERPRODUCT(u,y) index : (k) start : (s) ae abs(Ali, kK); imax i= i 


finish : (f); 
values, f; integer k,s,f; real u, v; 
‘comment INNERPRODUCT forms the sum of u(k) X 
‘ vik) fork=s,s+1,...,f.Ifs>f, the value 
of INNERPRODUCT is zero. The substitution 
of a very accurate inner product procedure 4: 


end 2; 

pivot[k) := imax; 

comment We have found that Alimax, k] is the largest 
pivot in column k. Now we interchange rows k and imax; 

ifimax # k then 

begin for j := 1 step 1 until n do 


i would make CROUT more accurate; 5. begi 

“begin _ rane . Pi 7 F 

i LEMP := Alkj};  Alk, j] := Alfimax,j]; 

; weal bh; Afimax,j] := TEMP 

oe 0; fork oe Ts h:=h+uXv; end 5: — 

; INNERPRODUCT := h APMHO ae: Glin de SE ARAM eels Ties eee 

sen INNERPRODUCT; ae := bik]; blk] := blimax];  blimax] := TEMTL 

“procedure CROUT (A, b, n, y; pivot, INNERPRODUCT); comment The row interchange is done. We proceed to the 

value n; array A, b, y, pivot; integer n, pivot; elimination; 
real procedure IN NERP RODUCT; if A{k, k] = 0 then go to singular; 

‘comment This is Crout’s method with row interchanges, as for i := k+l step | until n do 

" formulated in reference [Ll], for solving Ay = b begin quot := 1.0/A[k, k]; | Afi, k] := quot X Afi, k] 
and transforming the augmented matrix [A b] end; 
into its triangular decomposition LU with all for | := k+1 step 1 until n do 
L(k, k] = 1. Tf Ais singular we exit lo ‘singular,’ Alk, j] := Atkjl-— INNERPRODUCT (Aik, pl, 
a non-local label. pivot{k] becomes the current Alp,j], p, 1, k—1); 

Bo. row index of the pivot element in the k-th b[k] := bik] - INNERPRODUCT (AlLk,p], bipl, p, 

‘i column. Thus enough information is preserved 1, k-1) 

+ for the procedure SOLVE to process a new end 1; 
right-hand side without repeating CROUT. comment ‘The triangular decomposition is now finished, 
The accuracy obtainable from CROUT would and we do the back substitution; 
be much increased by calling CROUT with a for k := nstep —1 until I da 
more accurate inner product procedure than yik] :-= (b{k] — INNERPRODUCT(AIk,p], ylpl, p, 
INNERPRODUCT; k+1, n)/A[k, k] 

pein. end CROUT; 


integer k, i,j, imax, p; real TEMP, quot; 
fork := 1 step t until n do 

“1: begin 

‘TEMP := 0; 

fori := k step 1 until ndo 


procedure SOLVE (B, ¢, n, z, pivot, INNERPRODUCT); 
value n; array B, ¢, z, pivot; integer n, pivot; 
real proceedure INNERPRODUCT; 

comment SOLVE assumes that a matrix A has already been 
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transformed into B by CROUT. but that a new 
column ¢ has not been processed. SOLVE solves the 
system Ay = c, and the output « of SOLVE is pre 
cisely the same as the output y of the procedure 
statement CROUT (A, ¢. n, y, pivot, INNER- 
PRODUCT). However, SOLVE is faster, because 
it does not repeat the triangularization of A; 


begin 
integer k; real TEMP; 
for k := I step 1 until n do 
begin 
TEMP := clpivot(k]]; elpivat(k]] := elkl; cfk] c= 
TEMP; ef{k] := elk] INNERPRODUCT(BIK. pi, 
cip}, pl, k - 1) 
end; 
for k := nstep —l until 1 do 
2(k) c= (efk] — INNERPRODUCT(Bik,p], zfpl, p. 
k+1, n)/B[k, ki 
end SOLVE 
REFERENCE 
(1) J. H. Winkinson, theory and practice in linear systems, pp 
43-100 of Joun W. Carr IIT (editor), Application of Advanced 
Nuaerical Analysis to Digital Compuiers, (Lectures given at 
the University of Michigan, Summer 1958, College of En- 


gineering, Engineering Summer Conferences, Ann Arbor, 
Michigan. [1959]). 


ALGORITHM 17 

TRDIAG 

C. F. Spraeun HI 

General Atomic Division of General Dynamics Corp., 
San Diego, California 


procedure trdiag (a,b,c,d) order ; (a) result : (x); 
value n; array a,b, ¢, dx; integer n; 
comment this procedure gives the solution to the tri- 


diagonal system of linear algebraie equations: 
axe + bixi + di = 
aiXign + bixit CjX jin + dj —_ 0, i= 2.3, Pot ey n-1 
baXn + @nXaat + dn = 
‘This method is offen used to obtain solutions to 
second order difference equations; 


begin array gamma {1 :n—1]; imtegeri; real y; 
gamma [1] := ~—a{l]/b{l]; 
xf{lj := —d{t]/b{L]; 
for i:= 2step | until n—1 do 
begin y = bli] + efi] X gamma [i — 1]; 
gamma fi} = —alij/y; xfil:= —(¢fi} x xfi-1} 


+ dliji/y end; 
xin} i= —(efn] X x[n—-1] + d{n))/(b[n] + einf 
* gamma (n—1}); 
for i:=n step —| until 2 do 
xi — 1] := x{i] X gamma {i ~ 1] + xli ~ O] 
end irdiag 
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Covtributions to this department aust be in the form 
: stated in i Algorithms Department poliey 
£ (Communications, Vobruary, 


statement 
ope that ALGOL 69 
ce Conununications, May, 1960), 
Contributions should be seat tod. TW. Werstein, Computa. 


notation should be used 


tion Laboratory, Natioual Bureau of Standards, Washing. 
ton 25, D.C. Algovilhous should be in the Publication form 
of ALGOL 60 and written ina style patterned after the most 
recent algorithms appearing to this department, 
Although cach algorithm has heen tested by its gon 
iributor, no warranty, express or implied, is made by the 
contributor, the editor, or the Association for Computing 
Machinery as to the accuracy and functioning of the al- 
gorithm and related algorithni material and no responsi. 
bility is assumed by the contributor, the editor, or the 
Association for Computing Machinery in connection there- 


with. 

The reproduction of algorithms appearing in this de- 
partment is explicitiv permitted without any charge. When 
reproduction is for publication purposes, reference must be 
made to the algorithm author and to the Communications 
issue bearing the algorithm. 


ALGORITIIM 18 

RATIONAL INTERPOLATION BY CONTINUE 
FRACTIONS 

R. W. FLoyp 

Armour Researeh Foundation, Chicago, [linois 


comment This procedure fits to m given points (xi, yi) a col 
tinued fraction in the form 
art (xX) / (at (XX 9) / (aah (Xm Xa) / oe + (Sm Seat) /Am))+*)) 
Tt also simplifies the continued fraction to a rational functic 
(NoP Nix 4 Noes 8) (Dp Dix + Daogs 8), 
where deg is at most m + 
procedure confr(m,x,y,a,N ‘D); 
real array x,y,a,N,D; inveser m; 
begin real aa,xx,T; integeri,j,k; realarray P,Q|0O:m+ 
switch sw := swl, sw2; 
for j := lstep] until mdo 
begin aa := y{j]; xx := xij]; 


fori istep 1 until }—1 do 
aa :=(xx—xfip/(aa—ali}); alj] c= aa 
end; 
k:= 1; PO} := 1; Q{O) := afl]; 
mult :forj := lstep | until m + 2 do P{j| := Qij) = 9: 
for i := 2 step 1 until m do 


begin for j := i + 2 step —1 until | do 
begin T := ali] X Qf] — xfi— 1] X Pl) + PU 
Pili] QU Qlil = 
end; T := afi] X QQ) — xli-1] x PiO); 
PIO] := QO]; QO) := 
end; go to sw{k]; 
swl : for j := Getep Luntilm + ; do N{jl = QU) 
k= 2; Pl0}-= 0; Q0):= 1; go to mull; 
sw2 for} c= Datep 1 Luntilm + £246 Di} s= QU) 
end procedure 


